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The Tully-Fisher Relation for 25,000 SDSS Galaxies as Function of 
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ABSTRACT 

We construct Tully-Fisher relationships (TFRs) in the u, g, r, i and z bands and stellar mass 
TFRs (smTFRs) for a sample of 25, 698 late spiral type galaxies (with 0.045 < z < 0.085) 
from the Sloan Digital Sky Survey (SDSS) and study the effects of environment on the re- 
lation. We use SDSS-measured Balmer emission line widths, i>fwhm, as a proxy for disc 
circular velocity, v c hc- A priori it is not clear whether we can construct accurate TFRs given 
the small 3" diameter of the fibres used for SDSS spectroscopic measurements. However, 
we show by modelling the Ha emission profile as observed through a 3" aperture that for 
galaxies at appropriate redshifts ( z > 0.045) the fibres sample enough of the disc to obtain a 
linear relationship between wfwhm and i; c irc, allowing us to obtain a TFR and to investigate 
dependence on other variables. We also develop a methodology for distinguishing between 
astrophysical and sample bias in the fibre TFR trends. We observe the well-known steepen- 
ing of the TFR in redder bands in our sample. We divide the sample of galaxies into four 
equal groups using projected neighbour density (£) quartiles and find no significant depen- 
dence on environment, extending previous work to a wider range of environments and a much 
larger sample. Having demonstrated that we can construct SDSS-based TFRs is very useful 
for future applications because of the large sample size available. 
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1 INTRODUCTION 

The TFR is an observed correlation between the luminosity and 
disc c ircular velocity, t> c irc of disc type galaxies JTullv & Fisher! 
Il977t) . This fundamental empirical relationship reflects impo rtant 
physics in galaxy formation (e.g. iMo. Mao & White] 1 19981) and 
serves as a distance i ndic ator to galaxies (e.g. | Tulfy & FisheJ 19771 
iTullv & Pierc e 2000 and Freedman et al. 200l]). The luminosity of 
a galaxy is proportional to the stellar mass, while the rotation curve 
(RC) circular velocity is determined primarily by the dark matter 
halo (which is significantly more massive), so the TFR is essen- 
tially a relation between the dark matter component and the lumi- 
nous baryonic component of a ga laxy. 

The TFR was discovered bv lTullv & Fisher] d 19771) who mea- 
sured luminosity versus 21 cm neutral hydrogen emission line 
width with the 91 m National Radio Astronomy Observatory 
(NRAO) for a sample of 10 spiral galaxies. Using 21 cm data is 
advantageous because of the ability of a radio telescope to sam- 
ple the neutral hydrogen to very large galactic radii, obtaining 
data up to the plateau of the RC. Since its discovery, the TFR 
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has been investigated using other measures of disc circular ve- 
locity as wel l, such as optical Ha emission with long-slit spec- 
tros copy (e.g. Courteau 1997, Kannappan, Fabricant & Franxl2002l 



and lPizagno et alj|2007 ). TFRs have also been reconstructed from 
optical one-dimensional spectra (integrated line-of-s ight velocity 
widths) for galaxies up to z ~ 1 (Weiner et al. 2006b. a). 

There have been many studies of the TFR, each with its own 
selection criteria, but generally falling in two broad camps. The first 
camp is thos e interested in the phy s ical processes which give rise to 
the TF R , e.g.lPizagno et alj J2007h ;lK~annappan. Fabricant & Franx 



d2002h ; iBell & de Jond feOOlh : IPuechet alJ fcoioh 



Courteau 



( 1997) and many others. The introduction to ICourteau|[l997[ gives 
an excellent revie w. These surveys tend to have broader selec- 
tions, for example, ICourteaul dl997l) select bright (mj < 15.5), 
large (diameter > 4'), inclined (55° < i < 75°) galaxies 
from the UGC, eliminating only galaxies with large dust extinc - 
tion or peculiar/interacting morphologies dCourteau et al J [T993I) . 
Scatter in the TFR remains a large uncertainty in its use as a 
distance i ndicator, so physical explanations of this scatter such 
as that of iKannappan. Fabricant & Franxl d2002l) are a key area 
of research. The second camp is intereste d in the TFR a s a em - 
pirical law for measuring distances, e.g. Tully & Fisher ( 
Mathewson, Ford & Buchhom ( 1992); Mathewson & Ford 



1977); 



1996); 
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Springob et al. (2007) and many o thers. The use of t he TF R as 
a distance indi c ator i s reviewed in IStrauss & Willicki d 1995b and 
iTullv & Piercel d2000h . To minimise scatter in distance measure- 
ments, these studies tend to have narrower selection criteria, pref- 
erentially selecting late-type spirals (Sb - Sd) and luminosity pass- 
bands where intrinsic scatter is minimised. Despite differences, 
these two camps are closely linked: a better physical understanding 
of the TFR provides a better calibration of disk galaxy luminosities 
and hence more accurate distances. 

We will make extensive use of o ne such study here f or cal- 
ibration and comparison: the paper by iPizagno et all d2007l) . That 
work investigates a broadly selected input sample of 234 galax- 
ies drawn from the SDSS with —22 < M r < —18.5 to study 
the intrinsic scatter in the TFR. They obtain reliable Ha RCs from 
long-slit spectroscopic measurements with the Calar Alto 3.5 m 
telescope and the MDM 2.4 m telescope for 162 galaxies from 
the sample. This broadly selected sample does, however, contain 
typical magnitudes of M r < —18.5 and r-band scale heights of 
8", meaning it only samples the brightest and largest low-z galax- 
ies in the SDSS. They find the slope of the TFR systematically 
steepens from —5.5 ± 0.2 mag (log 10 km s -1 ) -1 in the g band 
to —6.6 ± 0.2 mag (log 10 km s -1 ) -1 in the z band. The intrin- 
sic scatter of the TFR is 0.4 mag in the g, r, i and z bands, and 
is attributed to being driven largely by the variations in the ratio 
of the dark to luminous matter within the disc galaxy population 
because correlations of the TFR with galaxy properties such as 
colour, size and morp hology were found to be weak. The slope of 
the TFR derived by IPizagno et alj is shallower than contemporary 
studies, which they describe as an effect of Malmquist-type biases 
for which they include no correction. We have chosen this sample 
for comparison here because it is based entirely on SDSS data. 

The large intrinsic scatter in the TFR suggests that the TFR 
may depend on properties external to the galaxy, such as the envi- 
ronment (i.e., the local number density of galaxies). Environment 
is known to have a strong impact on galaxy evolution, and, in the- 
ory, could have an effect on the TFR. For example, if galaxies in 
a higher density environment undergo accelerated star formation 
due to tidal interactions or mergers, then they may be overly lumi- 
nous for their rotational velocity and the TFR will be altered. Or, a 
galaxy fa lling into a cluster can have a s ignificant amount of its gas 
stripped dAbadi. Moore & Boweilll999l) . which can alter the TFR. 
Studying the TFR as a funct ion of environment h ence has impli- 
cations for galaxy evolution. iMasters et alj d2006l) conducted two 
simple tests for an effect of environment on the TFR and found no 
variation across a range of cluster environments in a sample of 5000 
galaxies as a function of clusters. 

The goal of this paper is the measure the TFR for a large sam- 
ple (~ 25, 000) of Sloan Digital Sky Survey (SDSS; lYork et alj 
2000) galaxies as a function of environment with very little sam- 
ple selection criteria except for a colour cut, which will se- 
lect mostly late disc type galaxies from SDSS. The TFR has 
not been investigated as a function of environment statistically 
before, although comparisons have been made between cluster 
and field TFRs as well as the TFRs for differen t clusters (see 
Pierce & Tullvll 19921 Iffliiz Rieke & Caldwellll20q3l. IMasters et alj 
20061 Springob et al. 2007, Springob et alj 120091 and references 
therein). iHinz. Rieke & Caldwelll d2003l) gather data on 13 Coma 
Cluster SO galaxies and 8 Virgo Cluster SO galaxies to compare 
their TFRs with the TFRs for late-type spirals and field SOs. They 
find no difference between the TFRs for field and cluster SOs, and 
only a small offset fro m the TFR for late-ty pe spirals. However, 
the sample size is small. [Masters et al] d2006t) studies ~ 5000 field 



and cluster spiral galaxies that have either 21 cm or Ha line width 
measurements to investigate the i band TFR. The sample comes 
from several datasets from the 1990s as well as significant amounts 
of new data, which ar e carefully homogenis ed, and is presented i n 
Springob et al d2007l) (see also the erratum Ispringob et al1l2009h . 



Masters et al J 1 2006) find a morphological dependence on the TFR: 



earlier type spirals are observed to be on average dimmer than later 
type spirals at a fixed velocity width, which results in a shallow- 
ing of the TFR slop e for earlier type spirals. This is attributed by 
Mas ters etaH d~2006) to be the reason why field versus cluster TFRs 
may be different (since early-type spir als are more l i kely to be 
found in over-dense regions). However, IMasters et alj d2006l) find 
no difference in the slopes, zero points, and scatters of the TFRs 
for 31 nearby clusters and groups, and see no trend in the TFR 
residuals with a galaxy's projected distance from their cluster cen- 
tre when investigating the impact of environment. 

In the present work, we directly investigate the TFR as a func- 
tion of the local density E, determined from the projected distances 
to the fourth and fifth nearest neighbours, ranging from voids to 
values typical for clusters for a broad sample of SD SS galaxies. 
The p rojected densities have been calculated as in iBaldrv et alj 
(2006), which have been applied in previous studies to exam- 
ine the environmental dependence of a num ber of galaxy prop- 
erties for large samples of SDSS galaxies. IBaldrv et alj d2006l) 
use E to find that the fraction of r ed galaxies increases with in- 
creasin g E as well as galaxy mass. IMouhcine. Baldrv & Bamford 
d2007l) investigate the stellar mass-metallicity (gas-phase oxy- 
gen abundance) relation in star-forming galaxies as a function 
of projected neighbour density and find no dependence of envi- 
ronment on the relation, and conclude that the evolution of the 
galaxies is largely independent of t heir environments. However, 
IMouhcine. Baldrv & Bamfor j d2007l) do find marginal increase in 
the chemical enrichment level at a fixed stellar mass in denser en- 
vironments, suggesting that environments can play a moderate role 
in galaxy evolution. 

The organization of the paper is as follows. In Section [2] we 
show by modelling aperture-biased Ha line profiles why we can 
expect to obtain a TFR using SDSS-observed line widths. In Sec- 
tion [3] we discuss the data and sample selection used to construct 
our TFRs. In Section [4] we investigate the environmental effects 
on the TFRs in our sample (and find no strong dependence on 
en vironment). We find our TFRs are consistent with the TFRs 
of IPizagno et alj ( 120071) . In Section[5]we discuss the implications 
of our results and our conclusions. The model details are found 
in the Appendix, as well as discussion of effects which may af- 
fect the sensitivity of our method (seeing, Malmquist bias, rota- 
tion curve shape). We use standard cosmological parameters of 
H = 71 km s" 1 Mpc" 1 , Q m = 0.27 and fi A = 0.73 throughout 
this paper. All errors are quoted at the 95 per cent confidence level 
unless otherwise stated. 



2 WHY WE CAN EXPECT TO OBTAIN A TFR WITH 
SDSS LINE WIDTHS 

We investigate a large sample of galaxies from the SDSS 
dYork et al. 2000). We use as a measure of the circular disc ve- 
locity the SDSS-observed, inclination-corrected Balmer emission 
full-width half-maximum (FWHM) velocity width, vfwhm (which 
have been reprocessed in the Max Planck for Astronomy and Johns 
Hopkins University (MPA-JHU) spectroscopic re-analysis). These 
widths suffer from aperture bias since SDSS galaxy spectra are ob- 
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tained through a 3" diameter fibre, which means that the fibres 
do not cover enough of the galaxy all the way out to the maxi- 
mal disc circular velocity v c i IC to provide direct measurements of 
n c i rc . However, we may still be able to obtain a TFR if v c i IC is 
correlated with the observed, aperture-biased ufwhm- A number 
of factors may aid in the existence of a relationship between i> c i rc 
and WfWHM. More luminous galaxies are farther away and hence 
smaller in angular size so the fibres cover more of the galaxy. Less 
luminous galaxies are nearby but physically smaller (in addition, by 
the TFR, the value for u c i rc is also smaller which will also decrease 
the discrepancy between the actual v c i IC and its value implied by 
bfwem)- Seeing (typical FWHM is ~ 1.4") may potentially help 
by increasing the observed FWHM, bringing it closer to the actual 
FWHM as would be observed by a large enough fibre, since some 
light entering the fibre from the galaxy comes from beyond the 
3" aperture projected onto the galaxy (although more light, from 
within the 3" aperture projection, is scattered out). 

Here, we develop a model for the observed Ha profile of a 
disc galaxy as measured through a 3" aperture, accounting for 1.4" 
FWHM seeing. The model assumes an infinitely thin inclined disc 
with an exponential surface brightness profile, a radial velocity pro- 
file characterised by an asymptotic circular velocity and a turn-over 
radius, and constant velocity dispersion across the galaxy. We then 
use our model to create a synthetic population with model parame- 
ters chosen according to distributions inferred from the SDSS and 
the local TFR. We investigate the relation between the true v c irc 
and the observed inclination-corrected wfwhm as a function of red- 
shift to learn at which redshifts it is possible to construct reliable 
TFRs using SDSS line widths. The details of the model are found 
in Appendix [A] An example Ha emission profile with and without 
the SDSS 3" aperture bias is presented in FigureQ] The modelling 
performed here in Section |2]is used to inform our final sample se- 
lection (Section [3) of disc galaxies from SDSS Data Release Eight 
(DR 8) l lAiharaet alj|201lh for studying the TFR. 

Previous studies in the literature have performed similar sim- 
ulations of line widths as observed by fibr e spectra t o unde rstand 
their observed distribution. For example, iRix et al | Jl997h sim- 
ulated the expected distribution of (observable) [Oil] emission 
line widths in galaxies (for a much smaller sample than ours) as 
would be seen by the AUTOFIB fibre spectrograph on the Anglo- 
Australian Telescope (AAT). They assumed various galaxy circu- 
lar rotation speeds, and included effects of random viewing angles, 
clumpy line emission, finite fibre aperture and internal dust extinc- 
tion on the profile of the emission line. 



2.1 Relation between v c i IC and dfwhm 

The free parameters in our model for Ha emission are inclination i, 
exponential scale radius h, dynamic turn-over radius r t and asymp- 
totic disc velocity u c i rc . Assuming distributions for these parame- 
ters that are typical for SDSS disc galaxies allows us to test for 
a correlation between v c i IC and ufwhm- A reasonable correlation 
would allow for the construction of a TFR using SDSS-measured 
line widths. The values for ufwhm we obtain are from the FWHM 
of the model profiles and inclination corrected (i.e., divided by 
sin(i)). Inclination-correction reduces the scatter between w c i rc and 

fFWHM- 

We choose inclinations i between 55° and 90°. Assuming ran- 
dom orientation, this is achieved by drawing a value for sin(i) from 
a uniform distribution between sin(55°) and sin(90°). The reason 
for this choice of inclinations is that in our sample selection, we 
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Figure 1. Example Ha emission profile calculated from our disc model for 
parameters i = 80°, h = 3", rt = 2 and v c - lIC = 200 km s — 1 (model 
details are found in Appendix [Aj- The black solid line shows the full pro- 
file without aperture effects. The blue solid line shows how the same profile 
would appear if observed with a 3" aperture and 1.4" seeing, like in SDSS 
observations. A Gaussian is fit to the aperture-biased profile, shown as a 
blue dashed line, which nearly coincides with the actual aperture-biased 
profile. The full profile, on the other hand, may show a double-peaked char- 
acteristic (depending on model parameters) as in this case. 



Table 1. Gaussian probability distributions for synthetic population param- 
eters 
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will want to only analyse galaxies that are more than 55° from face 
on, otherwise the projected line-of-sight velocities are too small. 

Establishing distributions for h, r t and v c i IC is more compli- 
cated. The distributions are a strong function of redshift because 
farther galaxies appear smaller, and more luminous galaxies can 
be seen farther away (and hence have higher v c i rc ). We establish 
the distribution from a sample of 44, 276 SDSS DR 8 galaxies, 
Sample A. The details of the selection of this Sample A are de- 
scribed in Section [3] Sample A consists of galaxies with redshifts 
0.005 < z < 0.085 and magnitudes -22.5 < M r < -15. 
We divide the sample into groups with redshifts centered around 
z = 0.01, 0.02, 0.08 and interval widths Az = 0.01. To ob- 
tain a distribution for h for a given redshift group, we assume that 
the Ha emission profile traces the profile of the emission of the 
band in which the Ha line is located. In the rest frame, Ha is at 
656.5 nm, which is in the r-band. Beyond z ~ 0.05, the emis- 
sion is redshifted to the i-band. We find that the exponential scale 
lengths have approximately a log-normal distribution, as shown in 
Figure [2] The probability distributions in each redshift bin are fit- 
ted with log-normal distributions, with parameters listed in TableQ] 
These distributions were used to construct our synthetic population. 

Next, we wish to describe a distribution for the velocity turn- 
over radius r t . This information is not measured by the SDSS. We 
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fitted with Gaussians, used 1 
using the local Tully Fisher relation found in lPizagno et alj d2007h : 
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to create our synthetic population. 



Mi = -6.321(log 10 wgo - 2.220) - 21.390 
= -6.321 log 10 v 80 - 7.357 
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(2) 



0.0 0.5 1.0 

Log 10 (h, I arcsec) 



Figure 3. Plot of rt versus h r for a sample of SDSS selected galaxies from 
Pizag no et all 120071) . 



investigate the samp le of 162 SDSS selected local galaxies from 
IPizagno et al.l J2007h with reliable Ha RCs obtained from long-slit 
spectroscopic measurements with the Calar Alto 3.5 m telescope 
and the MDM 2.4 m telescope. We find no significant correlation 
between r t and SDSS variables such as magnitude or h r (expo- 
nential scale radius in the r-band) which would help us predict a 
val ue of r t . A plot of r t versus h r is presented in Figure [3] In 
the IPizagno et al.l d2007» sample, the mean value of rt is 0.5 dex 
smaller than h r and has an approximately log-normal distribution. 
Therefore, given a distribution of h r , we assume a log-normal dis- 
tributions for r t with mean 0.5 dex smaller than the h r mean and 
the same standard deviation. 

Finally, to obtain distributions for i> c iro we use the distribu- 
tions of absolute i-band magnitudes M% and convert them to n c irc 



where t)go is the velocity (in km s 1 ) at a radius containing 80 
per ce nt of the i-band flux, with «go ~ « c i rc in the IPizagno et alj 
(l2007h sample. The distributions of v c i IC for the different redshift 
bins are plotted in Figure |4] which again are approximately log- 
normal, with best-fit parameters listed in Table Q] 

Using the above distributions for the free parameters, we con- 
struct a synthetic population of 300 galaxies for each of the red- 
shift bins using our disc model. We plot the relation between the 
actual v c i IC and SDSS-measured, inclination corrected dfwhm in 
Figure [5] We see that at low redshifts (z < 0.045), where galax- 
ies appear larger in angular size and hence we are probing only the 
central parts of the emission with the fibres, there is clear aperture- 
bias effects with the distribution of u c i rc versus ufwhm being steep 
and fan-shaped. At higher redshifts (z > 0.045), the relation can 
be modelled with a linear fit. At these redshifts, the residuals show 
no biasing pattern and the variance (sum of the squares of the dis- 
tance to the fit divided by the number of points) has dropped signifi- 
cantly: the variance from the low through the high redshift bins are: 
0.0147, 0.0126, 0.0086, 0.0073, 0.0055, 0.0045, 0.0041, 0.0035 
respectively. The scatter is still large in the relation due to aperture 
effects, however this does not prevent us from obtaining a good 
fit to the slope of the TFR (although this does make it difficult to 
use SDSS-measured line widths to study the intrinsic scatter in the 
TFR). 

The plots in Figure [5] can be thought of as indicating the ex- 
pected shape of the SDSS line width-based TFR (i.e., functions of 
wfwhm), since an assumed TFR (as function of u c irc) can be used 
to convert the j/-axis to absolute magnitudes. The shape of the dis- 
tribution is not expected to be linear if we select galaxies at all 
redshifts between 0.005 and 0.085, due to the aperture effects that 
become more prominent for lower z galaxies (which appear larger 
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Figure 4. Histograms showing the SDSS probability distributions of v c \ TC (derived assuming the local TFR from Pizagno et al. 2007), fitted with Gaussians, 
used to create our synthetic population. 



in angular size). The low z (and hence also dimmer) galaxies are 
expected to distinctly turn-off from the linear trend the higher z 
galaxies demonstrate due to the aperture bias. Therefore, in order 
to construct and analyse the TFR for a sample of SDSS galaxies, 
we will choose only galaxies with z > 0.045. 



We have chosen the lPizagno et al .1 feOOTt) relation because it is 
directly based on SDSS data, and therefore made for the simplest 
com parison. However, th ere are several caveats to this decision. 
The iPizagno et a 1. (2007) relation is shallower than most current 
TFRs, primarily because Malmquist-type biases are not adequately 
addressed (as described in their Section 5.2). There is also evi- 
dence for a morphological-type d ependence in the TF R, with later 
types having a steeper relation l lMasters et al.l l2006). By select- 
ing late typ e galaxies, th e true TFR in our sample may be steeper 
than in the Pizag no~et alj sample. Larger, but less b roadly selected, 
Malmquist corrected TFR samples are available (e.g lSpringob et all 
l200?tlMasters et alj2 006). The narrower selection of only late-type 
spirals in these larger samples matches better the colour selection 
we employ (see Section[3](. However, we do not include any addi- 
tional morphological pruning (e.g. for peculiar galaxies or mergers) 
common in these larger data sets. We also will not include a correc- 
tion for Malmquist-type biases in our sample (see Appendix |B3I >. 
Ultimately, our goal is to look for differences in the TFR as a func- 
tion of environment, so while the adoption of the Pizag nrJet alj 
TFR for our analysis may affect the overall TFR slope recovered, 
our method is self consistent and differences in slope with environ- 
ment will not be adversely affected by this choice of comparison 
TFR. We have checked that this is true by rerunning our a nalysi s 
using a steeper TFR (slope 7.85, based on [Masters et alj ((2006)) 
and find no difference in our conclusions (see Section 1431 . 



3 DATA AND SAMPLE SELECTION 

The data used for this investigation come from the SDSS DR 8. 
The SDSS is a photometric and spectroscopic survey mapping ap- 
proximately 7T steradians of the celestial sphere. Galaxies are im- 
aged in high quality in five bands (u, g, r, i and z), which pro- 
vide reliable magnitude measurements. SDSS galaxy spectra are 
obtained through a 3" diameter fibre. Initially, we select a large 
sample of 71, 933 galaxies with minimal amount of cuts using the 
online 'CasJobs' catalogue archive server SQL database interfac^B 
This sample, Sample 0, was selected according to the following cri- 
teria: 

(1) We selected galaxies using a colour cut with it — r < 2.22. 
The purpose of this colour cut is to filter ga laxies to choo s e ones 
with disc type morphology. The paper by IStrateva et alj (2001) 
finds that SDSS galaxies on a luminosity versus colour diagram 
have a bimodal distribution, with early type galaxies (redder) 
mostly having u — r > 2.22 and late typ e galaxies mostly hav- 
ing u — r < 2.22. Stra teva et alj ((2001) determine a dividing 
colour between the two populations at u — r — 2.22 by minimiz- 
ing the overlap of two Gaussians that define the two populations. 
This colour cut allows us to select mostly spirals, bu t we are ex- 
cludin g the red-end, dusty spirals as a consequence. Bald rv et alj 
d2004h show that dusty and bulge-dominated spirals may be found 
at u — r > 2.22 so these would be missing from our analysis. 
Therefore, our sample primarily consists of very late type (disc 
dominated) galaxies without significant bulges. Moreover, u — r 
colour is a function of inclination, with edge on galaxies appearing 
about 0.2 mag redder than face on galaxies, so this c olour cut cre- 
ates a slight preference for lower inclination galaxies (Mast ers et alj 
120101) . An alternative sample selection criteria one could use for 
future studies would be to use the visual classifications from the 
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Figure 5. The expected relation between i) c i rc and dpwhm assuming our simple disc model, free parameters drawn according to SDSS distributions, and the 
local TFR. The vfwhm are inclination-corrected. Since by the TFR 

w c irc corresponds to some magnitude Al (so the y-axes can be re-scaled as magnitudes), 
this plot shows how the shape of a scatter plot of absolute magnitude versus SDSS-measured velocity widths (ufwhm) looks like. Strong aperture bias is seen 
for galaxies at z < 0.045. The thick, red, dashed lines show the best-fit linear models to the correlation. The thin, black, dashed lines indicate below which no 
"fwhm' s should fall, given v c i IC and assuming an arctan velocity profile. 



Galaxy Zoo project l lLintottetal .1120081) . which could be used to 
selected a broad sample of disc type galaxies from very disc-like 
ones to very bulgy ones. However, our sample selection of very 
late type spirals does help disentangle environmental effects on the 
TFR related to bulges (bias could be introduced since there are 
more bulges in disc galaxies at hig her densities dDressletj|l980l 
iMouhcine. Baldrv & Bamfordl2007IV ). 

(2) We select galaxies with inclination above 55°. This incli- 
nation cut was made to choose only galaxies that are near edge-on 
rather than face-on in order to obtain reliable inclination-corrected 
velocity widths due to the rotation of the disc. At lower inclina- 
tions, the error on the line-of-sight velocity becomes comparable to 
the the velocity. The inclination of a galaxy was estimated from the 
measured axis ratio (b/a) for the galaxy using the formula: 

. ... Il - (b/a) 2 

an = \ ' , (3) 

v ; V 1-0.19 2 

which accounts for the finite thickness of the disc ( a discussion of 
this equation is found in lHavnes & Giovanellill 98*41 and the partic- 
ular fo rm of the equation we write here comes from lPizagno et alj 
d2007l) ). Axis ratios b/a are measured in five bands (ugriz). For 
each galaxy we average sin(i) as implied by the gri bands. We 
avoid using the u and z bands for inclination estimates since they 
are less sensitive to surface brightness features. We also do not just 
use a single band because a single band yields a histogram of the 
sin inclinations which shows a few narrow troughs and peaks on top 
of an otherwise flat distribution most likely due to the fitting algo- 
rithm, and averaging results gives a flattened histogram (as would 
be expected due to random inclinations). We select galaxies with 
sin(i) > 0.6 (corresponding to > 55° inclination) for our sample. 

(3) We select galaxies that have emission line measure- 



ments from the MPA-JHU spectroscopic re-analysis. Velocity 
widths (standard deviations) are mea sured simultaneously in all 
of the Balmer lines, as described in iTremonti et al. (2004) and 
Brinchmann et al. I d2004l) , where the full details may be found. 
These line width measurements are made with a more improved 
accuracy in continuum subtraction than the default SDSS spectro- 
scopic pipelines. 

(4) We only select galaxies that have good photometry with 
no important warning flags and accurately measured spectroscopic 
redshifts with no warning flags. 

(5) We make a redshift cut to select galaxies between 0.005 < 
z < 0.085. The reason for this is that only galaxies in this redshift 
range can have accurate projected neighbour densities (£) mea- 
sured by SDSS. A discussion of £ will follow shortly in this sec- 
tion. 

The data retrieved from SDSS for each galaxy are the fol- 
lowing: spectroscopic redshift z; model magnitudes in the five 
ba nds (ugriz) corrected for Galactic ex tinction computed follow- 
ing lSchlegel. Finkbeiner & Davisl ( n~998l) : axis ratios b/a in the five 
bands; K COII values in the five bands (from the Photoz table in 
CasJobs) used for cosmological correction in the calculation of ab- 
solute magnitudes; the standard deviation, <r raw , of the Balmer line 
widths; and total stellar mass es timates of the ga l axies (median 
estimates were used) derived in iKauffmann et alj d2003h. w here 
full details may be found (see als o iBrinchmann et alj |2004| and 
iBaldrv. Glazebrook & Drivei]|2008l) . photometry of the galaxies. 

Absolute magnitudes, M, in each band are calculated from the 
apparent magnitude m and redshift z according to 

M = m-51og 10 (£) L /pc)-5-iCcorr (4) 
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where Di, is the luminosity distance, which we approximate as: 
D L {z) = (1.831og 10 (l + z) + 0.375 sin(0.67z)) 
X (1 + z) X 4.23 x 10 6 kpc. 



(5) 



The raw line widths er raw (standard deviations) are converted to 
inclination-corrected FWHMs according to: 



«fwhm = 2.3548- 



(6) 



We ultimately wish to study the TFR as a function of envi- 
ronment, so we need a meas ure of environmental density around 
galaxies. iBaldrv et alj d2006l) calculate projected neighbourhood 
densities, S , for S PSS DR 4 galaxies. To estimate the density, 
Baldr v et al ] J2006I) determine for each galaxy the TVth nearest 
neighbour projected densities: 



N 



(7) 



where djv is the projected comoving distance to the iVth nearest 
neighbour to a galaxy that is a member of the density defining pop- 
ulation (DDP) of galaxies within ±Azc = 1000 km s" 1 . The DDP 
consists of galaxies with M r < —20 — 1.6(z — 0.05), that is, it 
consists of a volume-limited sample (for z < 0.085). Full details 
for the calculation of Ejv, including co rrections for galaxie s near 
the photometric edge, may be found in IBaldrv et al.1 ( 120061) . The 
be st estimate for the p rojected neighbourhood densities E is found 
bv lBaldrvetalJteOOeh to be: 



log 



ic,S = - (log 



l E 4 +log 10 E 5 ) 



(8) 



Baldry recalculates E for SDSS DR 7 galaxies, which are avail- 
able on his research pag^El We cross-match the spectroscopic IDs 
of these sources with our Sample from DR 8. We obtain a new 
sample, Sample A, of 44, 276 galaxies with density measurements, 
keeping only galaxies with ufwhm > 50 km s" 1 . Sample A was 
the sample we used in the analysis of Section [2] in order to inves- 
tigate at which redshifts might we be able to obtain a TFR with 
SDSS line widths. 

From Section|2]we have learned that it is best to use a sample 
of galaxies with z > 0.045 to reduce aperture biases and con- 
struct an SDSS-based TFR. The study bv lKewlev. Jansen & Gellej 
d2005h have also found that using the SDSS spectra of the inner 
parts of galaxies at z < 0.04 can have strong systematic effects 
in the estimation of their global properties such as star formation 
rate, metallicity and reddening. We create Sample B, consisting of 
25, 698 galaxies, by eliminating galaxies with z < 0.045 from 
Sample A. Sample B is the final sample we construct, which we 
use for studying the TFR as a function of environment. Galaxies in 
Sample B have typical magnitudes —22.5 < M r < —18.5. 



4 RESULTS 

The aperture-biased emission profile modelling we performed in 
Section[2]explains why we can expect to obtain a TFR using SDSS 
line widths measured by a 3" fibre. In this section, we demonstrate 
that we do obtain linear TFRs using the actual data, which makes it 
possible to study the TFR as a function of environment. 

Firstly, we plot absolute magnitude in the i-band, Mi, versus 
ufwhm for Sample 0, presented in Figure [6] This sample contains 



2 http://www.astro.ljmu.ac.uk/~ikb/research/ 
© 201 1 RAS, MNRAS 000. mmi 



-16 
-14 




Log 10 (vfwhm / km s ) 

Figure 6. Plot of absolute i-band magnitude versus inclination-corrected 
dfwhm for a sample of 71, 933 SDSS galaxies. Galaxies with z < 0.045 
are in red, and higher redshift galaxies in blue. Over-saturated points are 
plotted as log-scale contours. This is a Tully-Fisher plot, except that there is 
a turn-off due to strong aperture bias at low redshifts z < 0.045 (red). The 
distribution shows the same characteristic distribution as predicted by our 
disc emission models in FigurefS] For comparison, the pl ot also shows the i- 
band TF R as a function of vgp f rom lPizagno et al.1 1 2007) (long-dashed line) 
and from Masters et al.1 (2006) (short-dashed line). The average z for the 
data clump of galaxies that turn off from a linear trend at M; > — 17 mag 
is z = 0.011. 



galaxies with redshifts 0.005 < z < 0.085, so we expect to see 
strong aperture effects from the low-redshift galaxies. The purpose 
of Figure [6] is to see where the galaxies lie and compare its gen- 
eral shape with the plots of Figure [5] which shows the predicted 
shape of the TFR from our simple disc models. We do see the low- 
z (hence low luminosity) turn-off of galaxies, and we are justified 
in excluding galaxies with z < 0.045 in Sample B. There is also 
some degree of turn-off expected for low-lu minosity galaxies wit h 
«circ < 90 km s _1 due to physical reasons (McGaugh etal. 2000). 
This turn-off is due to the fact that these low luminosity galaxies are 
very rich in gas, and if the gas mass in the disc is added to the stellar 
mass, then the turn-off can be rectified and these galaxies may be 
returned to the TFR trend line. However, we excluded these low- 
luminosity galaxies in Sample B when we made our redshift cut, 
and therefore do not have to worry about correcting for this effect. 

We plot the SDSS-based TFRs in five bands (ugriz) for 
Sample B in Figure |7j We obtain linear trends with large scat- 
ter (~ 0.5 mag) due to the aperture effects as predicted in Sec- 
tion^ The estimated additional scatter by aperture bias is ~ 0.3- 
0. 4 mag, comparable to the 0.4 mag intrinsic TFR scatter found 
in IPizagno et al.1 (20071) . Linear models are fitted to data in the 
five bands, with equations listed in Table [3] We perform forward 
fits, and have not included errors on individual points in fitting 
the relation, an approach justified due to the uniformity of the 
SDSS dataset across such a large sample of galaxies (i.e. the er- 
rors on each point are effectively the same, so including them 
does not affect the result, but does greatly increase the process- 
ing complexity of the problem). We see steepening in the TFR 
from —2.53 ± 0.06 mag (log 10 km s -1 ) -1 in the u band to 
—3.87 ± 0.06 mag (log 10 km s -1 ) -1 in the z bands (it is ex 



pecte d for the TFR to steepen in redder bands (Tulfy & Pierce 



| 200fj| : IVerheiienl I200TI; iKarachentsev et al.1 120021: 



Courtcau ct al 



20071 IPizagno et alj 120071 : iMasters, Springob & Huchral 12008)) 
There is no significant difference in the scatter between the differ- 
ent bands, since the aperture effects significantly contribute to the 
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scatter. In an ordinary TFR, one may expect l ess scatter in the red- 
der bands dAaronson. Huchra & Mould! 1979tlTullv & Piercel2000l : 
iMasters. Springob & H uchra 20081), where the mass to luminosity 
ratio of galaxies is more homogeneous since the redder bands are 
less sensitive to star formation and more sensitive to the older pop- 
ulation of star s which dominate the mass of the galaxy (although 
Pizag noet al.l d2007t) find the same intrinsic scatter in the g, r, i 
and z bands, possibly because of the fairly small sample size). In 
Figure|7]we also plot the smTFR for the sample, and the best linear 
fit to describe log 10 (A / l* / 'Mq) versus log 10 (i)FWHM/km s _1 ) is 
found in Table [3] 

For simplicity we have not included an inclination depen- 
dent extinction correction in our analysis. Extinction due to dust 
in the disks of spiral galaxies is a strong function of viewing angle 
jMasters et ai]l2010l : Mailer et alj|2009r) . IMasters et all d2010h find 
extinction of a galaxy viewed edge on is 0.7, 0.6, 0.5, and 0.4 mag 
greater tha n for a face on view for the Sloan ugri passbands, re- 
spectively. Malle retal] show that the extinction is a linear function 
of log 10 b/a, where b/ais the observed axis ratio of the galaxy. For 
galaxies with 55° < i < 90°, this implies an extinction range in 
e.g. i-band of 0.13 to 0.4 mag. While this range is significant, it 
is still smaller than the overall scatter of ~ 2.5 magnitude. Future 
analysis of SDSS data in this way would benefit from including 
an inclination dependent correction, however, since the distribution 
of inclinations should not correlate with environment, the omission 
here does not affect our conclusions. 



4.1 TFR as function of environment 

To investigate the environmental dependence of the TFR, we di- 
vide Sample B into four equal-sized groups based on projected 
density quartiles. The minimum, maximum, and quartile values of 

log 10 (E/Mpc- 2 ) are: 

{min, Qi, Q 2 , Qs, max} = {-2.06, -0.74, -0.43, -0.06, 1.96}. 

(9) 

Call the four groups Gl, G2, G3 and G4, with Gl (see Table [2j 
corresponding to the galaxies with the least dense environments 
typical for galaxies in void-like regions and G4 corresponding to 
galaxies in the most dense environments typical to cluster galax- 
ies and galaxies at the peripheries of clusters. Here, we investigate 
whether there is a significant difference in the TFRs of galaxies 
belonging to these four different groups. 

We plot the SDSS fibre-based TFR as a function of wpwhm 
for each group in each of the bands in Figure [7] The best-fit linear 
models (forward fits) for the TFRs are listed in Table|3] In Figure|7] 
we also plot the SDSS-based smTFRs for each of the four density 
groups, with best fits listed in Table[3] The 95 per cent confidence 
(statistical) errors on the slopes and intercepts from 100 bootstrap 
realizations are included in the tables in parentheses (i.e., to obtain 
errors, for each realization the data set is re-sampled with replace- 
ment and fit with a linear relation to obtain a probability distribution 
for the slope and intercept). The slopes of these TFRs as a function 
of dfwhm are approximately half the value of the slopes of real 
TFRs as a function of v c i Ic . This is attributed to the large scatter 
and flattened slope in the relation between v c - 1TC and ufwhm due to 
aperture bias. We show shortly in this section with the aid of sim- 
ulations that these flattened fibre-based TFR slopes are expected if 
we assu me that the rel ation between magnitude and v c i IC are the 
TFRs of lPizagno etafl fe007l) . 

First, we notice that the TFRs for the four environmental den- 
sity groups are similar to each other for each of the five bands. 



Table 2. Definitions of density groups 



group 


log 10 (S/Mpc- 2 ) 


Gl 


[-2.06,-0.74] 


G2 


[-0.74,-0.43] 


G3 


[-0.43,-0.06] 


G4 


[-0.06,1.96] 




2.2 2.3 2.4 

Log 10 (vpwHM/kms" 1 ) 

Figure 8. TFR in the five bands for all SDSS galaxies in Sample B as well as 
the different density groups. The colours violet, green, red, black and blue 
correspond to it, g, r, i and z bands respectively. The thick dashed lines 
show the TFR for all the galaxies in Sample B. The shaded areas represent 
the 99.99 per cent mean prediction bands for the relationship. The thin lines 
show the TFR for the different density groups for the five bands. The TFRs 
for the least dense and most dense groups lie outside the shaded confidence 
region in all five bands. 



There does appear to be a small amount of steepening of the rela- 
tions with denser environments. The TFR steepens 11.8, 9.2, 6.5, 
5.2 and 4.3 per cent in the it, g, r, i and z bands respectively from 
the least dense to most dense groups. The relative steepening is 
larger in bluer bands. In Figure [8] we plot the TFR in each of the 
five bands for the entire sample, showing (statistical) 99.99 per cent 
mean prediction bands for the relation, and plot the TFRs derived 
for each density group. The TFRs for groups G2 and G3 lie within 
the confidence region, but the TFRs for groups Gl and G4 lie out- 
side. At this point, we devise several tests to investigate whether 
this steepening is significant and whether it is due to the environ- 
ment. 

To test the significance of the differences in slope, we ran- 
domly permute the densities of the galaxies and repeat dividing 
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Figure 7. TFRs in five bands and smTFRs for SDSS all galaxies (Sample B) and galaxies divided into groups by density. Over-saturated points are plotted as 
log-scale contours. Gl represents the group of galaxies with the lowest 25 per cent of environmental densities and G4 represents the highest 25 per cent. The 
slopes and scatter in each band appear mostly independent of environmental density, with slight increase of slope with higher density in the bluer bands. Thick, 
red, das hed lines show the be st-fit linear models. Thin black lines show the griz band relations expected from simulations (see Section l4~Tl and assuming the 
TFRs of lPizagno eta?] ilOOH) . 



Table 3. TFRs in five bands and smTFRs for the entire sample and environmental density groups (as function of ifwhm) 



u g r 



All 


-13.08(0.14) - 


2.53(0.06)i 


-13.31(0.13) - 


2.95(0.06)1; 


-12.88(0.13) 


- 3.36(0.06) 


Gl 


-13.30(0.29) - 


2.42(0. 13)u 


-13.51(0.23) - 


2.85(0.11)i; 


-13.01(0.25) 


- 3.28(0.12) 


G2 


-13.29(0.27) - 


2.45(0.12)u 


-13.49(0.23) - 


2.88(0.11)1- 


-13.04(0.24) 


- 3.29(0.11) 


G3 


-13.11(0.29) - 


2.53(0.13)1- 


-13.39(0.23) - 


2.93(0.11)1- 


-12.96(0.31) 


- 3.33(0.13) 


Gl 


-12.70(0.30) - 


2.71(0.13)i! 


-12.98(0.30) - 


3.11(0.14)1- 


-12.60(0.28) 


- 3.49(0.14) 


All 


i 

-12.56(0.14) - 


3.64(0.06)i 


2 

-12.26(0.14) - 


3.87(0.06)i; 


sm 

5.28(0.08) + 2.11(0.04)i 


Gl 


-12.66(0.26) - 


3.58(0.13)u 


-12.35(0.28) - 


3.81(0.13)1 


5.23(0.13) 4 


2.12(0.06)i 


G2 


-12.71(0.25) - 


3.58(0.11)i> 


-12.42(0.28) - 


3.80(0.13)i 


5.31(0.15) + 


2.09(0.07)1 


G3 


-12.65(0.32) - 


3.60(0.14)1- 


-12.36(0.32) - 


3.84(0.15)i 


5.33(0.15) + 


2.09(0.07)1 


G4 


-12.31(0.34) - 


3.76(0.15)u 


-12.05(0.23) - 


3.98(0.15)1 


5.31(0.18) 4 


2.10(0.08)i 



© 2011 RAS, MNRAS OOO.fTlfBI 



10 P.Moczetal. 



Table 4. TFRs for density groups assuming parallelism 



band 






TFR 






it 


-13.09 


- 2.53?) 


4- 0.0282:i 


- 0.02222 


- 0.01923 


9 


-13.34 - 


- 2.95ii 4 


0.04821 - 


- 0.005522 


- 0.008523 


r 


-12.91 - 


- 3.35i> 4 


0.067zi 4 


- O.OO8822 


+ O.OO2323 


i 


-12.59 


- 3.64d - 


h 0.07321 - 


f 0.012z 2 - 


f 0.005123 


z 


-12.30 


- 3.87ti - 


h 0.07721 - 


f O.OI822 - 


f 0.009623 


stellar 


5.32 4 


2.10d - 


0.052zi - 


0.024Z2 - 


0.01923 



the galaxies by density quartiles and fitting each group with a lin- 
ear model for 2500 realizations. For each realization, we compute 
the steepening ratio between the highest and lowest sloped groups. 
From the distribution of the ratios, we find that the steepening in the 
u band is significant at the 0.01 probability level and the steepening 
in the g band is significant at the 0.02 probability level. Steepening 
in the other bands are not significant at the 0.05 probability level. 
TFRs assuming that the lines are parallel in each of the four density 
groups are presented in Table |4] The variables Z\, z%, and 23 are in- 
dicator variables, such that Zi = 1 for a galaxy if it is in group Gi, 
else Zi = 0. In the u and g bands, where parallelism is rejected at 
the 3<r level, there is no coherent pattern in the y-intercept between 
the density groups assuming parallelism. In the r, i, and z bands 
the TFRs marginally brighten with increasing density, assuming the 
relations are parallel across the density groups. The smTFR shows 
marginal increase in stellar mass at fixed velocity width with in- 
creased density. 

A difference in the linear relations between magnitude and 
bfwhm for two groups of data may not necessarily imply that the 
TFR is not the same for the two groups. As we learned from our 
simulated population and FigurefJ] the slope is affected by the size 
distribution of the galaxies. Differences in size distributions in the 
two groups may therefore cause the difference in slope. Thus we 
create synthetic samples of galaxies (same sample size as the ob- 
served sample) for each of the four density groups, with model pa- 
rameters drawn from the observed distribution of sizes (and again 
we assume the local TFR to estimate a distribution for disc veloci- 
ties). The purpose of the simulations is to see whether they predict 
the same slopes for the relation between magnitude and wfwhm as 
observed for the four density groups, and also to obtain errors on 
the slopes. The distributions of parameters for galaxies in each of 
the four density groups are listed in Table [5] where it can be seen 
that there are some small differences which may affect the slope. 
The distribution of size s at a fixed stellar mas s is independent of 
environment in general ( Kauffma nn et al.l2004l) . however the slope 
of «circ versus dfwhm may be sensitive to small discrepancies in 
sizes. We fit the relation between n c i rc and bfwhm along with er- 
rors predicted from 30 synthetic sample realizations. The results are 
shown in Table[6] which indicates slight steepening in v c i rc versus 
dfwhm with density (and tight errors) and says that the observed 
SDSS-based TFR slopes should be around 1/2 the values of the 
real slopes. Next, we convert the v c j r c to griz magnitudes using 
the TFRs found by Pi zagno et alj 1 20071) and add random Gaussian 



scatter of 0.4 mag dPizagno et alj d2007l) find for their sample of 
162 galaxies that the TFR residuals are approximately Gaussian). 
We use ugo = u(3.03/i) = 0.80t> c h-c in the conversion of u c i rc to 
magnitude for our exponential-profile disc models. Doing so, we 
have created synthetic magnitude versus «fwhm plots just like the 
observed data and can use it to predict differences in the relation 
in the four density groups as well as better estimates of error in the 
linear fit parameters than with bootstrapping (see Figure |9] . The 
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Figure 9. The observed i-band data (blue) with the best fit SDSS-based TFR 
(thick, red dashed line) compared to simulated data (green; 800 data points 
are shown for clarity) and its b est-fi t linea r relation (thin, black line) which 
assumes the TFR o f lPizagno etalU2007l) and the galaxy size distribution 
of our sample. There is very good agreement between the two, despite the 
simplifying assumptions we have made. A histogram of the dfwhm for the 
real and simulated data is also shown, showing close agreement, but slightly 
more dispersion in the real data than the simplified model. 



SDSS-based TFRs for our synthetic population are presented in Ta- 
bleland plotted along with the real data in Figure |7]as black, solid 
lines in the griz bands. The predicted and observed slopes and in- 
tercepts agree remarkably well given the simplifying assumptions 
we have made. From the simulations we also obtain a good sense of 
what the errors on the parameters of the fits to the real data should 
be, and find that they are comparable to the bootstrap errors we pre- 
viously found. Notably, we predict steepening on the order of 7 ± 4 
per cent of the SDSS-based TFRs from the lowest to the highest 
density group in all the bands due to size distribution differences in 
the groups (see Figure [ITU. 

The agreement between our simulated data and observations 
suggests that there is no environmental effect on our TFRs. We can 
ask now how sensitive is our method, that is, how much would the 
slope of the TFR as a function of D c i rc have to change so that our 
synthetic populations would predict an SDSS-based TFR that lies 
outside the 95 per cent confidence interval for the observed SDSS- 
based TFRs. To ans wer this, we perturb the slopes of the TFRs of 
Pizag noet alJj2007T) . which we have used in our simulations, about 
log 10 (i) c i rc /km s _1 ) = 2.2 and repeat the simulations to predict 
the slopes of the observed SDSS-based TFRs. We find that we need 
to perturb the slopes by < 3 per cent so that the predicted fibre- 
based TFRs lie outside the 95 per cent confidence intervals of the 
observed relations in Table [3] Similarly, we find that we can per- 
turb the intercept of the TFR by < 2 per cent so that the predicted 
fibre-based TFRs lie outside the 95 per cent confidence intervals 
of the observed relations. Therefore, we can say that if there are 
environmental effects on the TFR in our sample, then it would be 
very subtle, having less than 3 per cent on TFR parameters. This 
3 per cent variation is comparable to the errors on the observed 
TFRs in the first place. 



4.1.1 Another test for environmental effects on the TFR across 
bands 

Another test we can perform to show that the steepening in the 
TFRs from the least to most dense groups is not due to environmen- 



© 2011 RAS, MNRAS 000,[T]|l5] 



TFR as function of environment 1 1 



a 

(X 



o 

II 

(X 



Figure 10. The simulated population of galaxies using the observed param- 
eter distribution predict a steepening of 7 ± 4 (red dashed line and shaded 
region) per cent in the slope of the SDSS fibre-based TFR from the lowest 
to the highest density group in all magnitude bands, due to small varia- 
tions in the size distributions of the different density groups (assuming no 
trends with density). The observed data agrees, within error, with this pre- 
diction. Without aperture biases, no steepening would be expected (thin, 
black dashed line) if there were no environmental trends with density. 



tal effects is to look at the observed ratio of steepening across dif- 
ferent bands and compare it across bands. An important realization 



is that given two groups of log^ 



versus log 10 wpwhm data 



which may have slightly different slopes si and S2, if the y-axis 
is converted to magnitudes using the real TFR in order to obtain 
magnitude versus log 10 vfwhm plots (such as the observed data), 
then the ratio of the new slopes s[ , s' 2 of the two groups should be 
the same in any magnitude band. This is because if the TFR is of 
the form M = a\og 10 (v) + b then: 



£i 

S2 



Ai(lo gl 



=)/Ai(log 

io ^fwhmJ 



A2(log 10 U c irc)/A2(log 10 HFWHm) 

aAi(log 10 u c irc)/Ai(log 10 dfwhm) 
aA 2 (log 10 u circ )/A2(log 10 wfwhm) 
AiM/Ai(log 10 dfwhm) 
A 2 M/A 2 (log 10 »fwhm) 



(10) 



(11) 



(12) 



(13) 



This gives the most powerful test to see whether the TFR alters 
with density because we can simply compare the ratios of steep- 
ening across bands and do not have to rely on the simple model, 
or a particular assumed "real" TFR, or trying to correct Dfwhm to 
SOme ZJcirc . The simulated data though is necessary to give us an 
estimate on the error/variation allowed on the slope ratios for the 
TFR to be parallel in all the density subgroups (but the magnitude 
of the errors is not expected to vary significantly with model pa- 
rameters). We have observed a slope increase in the actual data of 
11.8, 9.2, 6.5, 5.2 and 4.3 per cent in the ugriz bands, while the 
simulated data predict a steepening of 7 ± 4 per cent (see Table[8}. 
This is consistent with no dependence of the TFR on environment. 
A similar test for the intercept of the TFR finds no dependence on 
environment (see Table [8}. 

In summary, the simulations show that steepening is expected 
with environment in the SDSS-based TFRs due to the small differ- 
ences in size distributions for the different density groups, and most 
importantly that the observed slopes and int ercepts of the SDSS- 
based TFRs are consistent with the TFRs of |Pizagno~e t al. ( 2007) 
at the 95 per cent level (compare the red and black lines in Fig- 
ure [7}. 



Table 5. Gaussian probability distributions for observed parameters in den- 
sity groups 







arcsec) 


l°Slo( v c 


rc /kms 1 ) 


group 




(j 


A* 


a 


Gl 


0.521 


0.178 


2.08 


0.101 


G2 


0.524 


0.173 


2.10 


0.102 


G3 


0.522 


0.171 


2.10 


0.103 


G4 


0.521 


0.164 


2.10 


0.105 



Table 6. t; c i rc versus i>fwhm from simulations 



group 



logio^circ) 



All 0.573(0.008) log 10 (u PW HM) + 0.923(0.017) 

Gl 0.556(0.012) log 10 (u F WHivi) + 0.953(0.025) 

G2 0.563(0.017) log 10 (u F WHivi) + 0.946(0.035) 

G3 0.569(0.011) log 10 (u F WHivi) + 0.936(0.023) 

G4 0.595(0.015) log 10 (i> F WHM) + 0.88(0.031) 



4.2 Assuming a Different TFR in the Analysis 

We have already di scussed the limitations of the choice of the 
Pizag noet alj d2007l) TFR (Section l2.lt . Here we investigate the 
effects of using a different choice for the TFR for the analysis. 
We pick a steeper i-band TFR of Mi = —7.85 log ln vgp — 4.00. 
This T FR was obtained by increasing the slope of the Pizag ncTet alj 
j2007l) TFR to -7.85 about log 10 v so = 2.2. This steeper slope of 
7. 85 is based on the J- band TFR for late-type spiral (Sc) galaxies 
of Masters et al. 



i 2006h the correction between /-band and SDSS 



Jordi, Grebel & Ammonll2006l) ) 



i-band is small 

We repeat all of the procedures with this new TFR, includ- 
ing estimating the distribution of v c [ IC from the SDSS magni- 
tudes. This analysis is not equivalent to assuming a fixed inferred 
distribution of galaxy rotation profiles (hence v C j rc and ufwhm) 
and investigating the effect of altering the TFR on the observed 
SDSS DFWHM-based TFR (which would predict an increase in the 
slope of the observed SDSS-based). Re-inferring the distribution 
of Wcirc from the SDSS magnitudes using a steeper TFR leads 
to a narrower distribution of v c irc, which can flatten the slope of 
the correlation between w c i rc and dfwhm- With this steeper slope, 
the predicted SDSS-observed ufwhm -based TFRs for the den- 
sity groups Gl, G2, G3, G4 are: M t = -13.01 - 3.30v PW hm, 
M l = -13.03 - 3.31ufwhm, Mi = -12.82 - 3.42« PW hm, 
Mi = — 12.63 — 3.52i>fwhm respectively. The size of the errors in 
the slope and intercept are similar to the previous simulations. The 
predicted slopes are slightly shallower than observed in the SDSS 
data. However, the simulations still predict steepening of 1.07 for 
slope G4 / slope Gl. 

A steeper TFR does not necessarily mean that the predicted 
SDSS -ufwhm -based TFR will be steeper, if one has to infer the 



Table 8. Steepening due to differences in size and velocity distributions in 
density groups 





slope G4 / slope Gl 


intercept Gl / intercept G4 


predicted 


1.07 ±0.04 


1.03 ±0.02 


obs. u 


1.118 ±0.08 


1.047 ±0.03 


obs. g 


1.092 ±0.07 


1.041 ± 0.03 


obs. r 


1.065 ±0.06 


1.032 ± 0.03 


obs. i 


1.052 ±0.06 


1.027 ±0.04 


obs. z 


1.043 ±0.05 


1.024 ±0.03 
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Table 7. Expected TFRs (as function of dfwhm) from synthetic populations assuming simple emission model 



All 


9 

-13.06(0.09) - 


3.13(0.04)i; 


r 

-12.89(0.12) - 


3.41(0.06)i; 


i 

-12.60(0.14) - 


3.62(0.07)i; 


2 

-12.28(0.15) - 


3.77(0.08) 


Gl 


-13.19(0.21) - 


3.06(0.10)i; 


-13.01(0.22) - 


3.30(0.11)1; 


-12.74(0.22) - 


3.52(0.11)1; 


-12.47(0.23) - 


3.65(0.11) 


G2 


-13.15(0.22) - 


3.10(0.11)1; 


-13.01(0.27) - 


3.26(0.13)1; 


-12.74(0.32) - 


3.56(0.16)1; 


-12.41(0.31) - 


3.71(0.15)- 


G3 


-13.15(0.21) - 


3.10(0.10)i; 


-12.93(0.18) - 


3.40(0.09)i; 


-12.67(0.19) - 


3.59(0.09)i; 


-12.36(0.24) - 


3.74(0.11> 


G4 


-12.84(0.16) - 


3.25(0.07)i; 


-12.65(0.22) - 


3.54(0.11)i; 


-12.33(0.28) - 


3.76(0.13)i; 


-11.99(0.24) - 


3.92(0.12) 



distribution of i) c irc from the SDSS data. This is because the corre- 
lation between i; c i rc and ufwhm may flatten if there is a narrower 
distribution of fdrc- This slope may be changed/strengthened if 
we change the intercept of the assumed-TFR as well. Exploring 
the parameter space further, we find that a steep (slope 7.85) TFR 
of Mi = —7.85 log 10 i>8o — 4.25 is consistent with the slopes of 
the SDSS ufwhm -ba sed TFRs. Thi s TFR is equivalent to chang- 
ing the slope of the IPizagno et all d2007t) TFR to -7.85 about 
log 10 wso = 2.0. 

It appears that the right choice of two different slopes and in- 
tercepts may conspire in certain cases to give the same predicted 
SDSS «FWHM-based TFR. However, there is another prediction 
made by the simulations, namely the distribution of the observed 
vfwhm- A steeper TFR predicts a narrower distribution of v c i rc 
and consequently i>fwhm- So a useful check to make sure there 
is no conspiracy is to investigate whether there are significant dif- 
ferences in the distribution of the observed distribution of dfwhm 
across density groups (these differences we find, would be much 
larger than those caused by the small differences in size distribution 
across density groups). We find none: the mean and standard devia- 
tions of log 10 vfwhm are 2.21 and 0.12 in all density groups (with 
variation < 1 percent). Based on ou r simulations, a TFR slope of 
6.32 (based on Pizagno et aid (120071) ) and a steeper TFR slope of 
7.85 (based on lMasters et all ( 1200611 ) which predict the same SDSS 
iifwhm -based TFR, would predict differences in mean and stan- 
dard deviation in ufwhm at the level of 8 per cent. 



5 DISCUSSION AND CONCLUSIONS 

We have demonstrated that we can use SDSS-measured line widths, 
which suffer from 3" fibre aperture bias, to construct TFRs for 
galaxies at 2 > 0.045, where enough of the galaxy fits inside the 
fibre so that we can recover its global properties. Hence, the SDSS 
database can be very useful for studying the TFR on a large sam- 
ple of galaxies, and one of the best applications is to look at the 
effects of environmental density on the TFR. We did so, and found 
the following: 

• We divided our sample of galaxies into four groups based on 
density, and found slight steepening of the SDSS-based TFR slope 
with density in each of the five bands (ugriz), but this was at- 
tributed (using simulations) to be a consequence of aperture biases 
affecting the density groups, which had small variations in size dis- 
tributions. We conclude that there is no strong or statistically signif- 
icant environmental effect at the 95 per cent level on the TFR in our 
sa mple of 25, 698 gala xies and the data is consistent with the TFR 
o APizaeno et al\ ^200^) . The environments we investigated ranged 
from void-like regions, with typical £ of 0.05 Mpc -2 , to centres of 
galaxy clusters, with typical E of 20 Mpc~ 2 . Our method is sensi- 
tive enough to allow us to conclude that if there is an effect on the 
TFR due to environment, the change in the TFR slope or intercept 
is less than 3 per cent. 



• Steepening in the TFR is seen on the order of 4 to 11 per cent 
with increasing density, with the relative steepening being larger 
in the bluer bands. Our simulated populations of model galaxies 
suggest a steepening of 7 ± 4 per cent across all bands, account- 
ing for the size change in the groups. One concern, however, may 
be that perhaps some of the steepening is due to contamination by 
early type galaxies in our sample, which can have a bias ing ef- 
fect since they tend to be found in denser environments dDresslej 
1980, Mouhcine. Baldrv & B amford 2007) and we had only made 
a simple colour cut to select disc type morphology and may have 
included some of the less-red ellipticals in our sample. However, 
repeating the analysis with a stricter colour cut of u — r < 2.00, 
we witness the same degree of increase in the TFRs with increasing 
density (11.8, 9.1, 6.5, 5.2 and 4.3 per cent in the ugriz bands). 
In addition, there is no significant difference in the concentration 
indices Ci — Rgo / R50 of the galaxies in the four density groups 
(-R90 and i?5o are the radii enclosing 90 per cent and 50 per cent 
of the i-band Petrosian flux). The concentration index is a measure 
of morphology - the significance of the bulge - for galaxies, with 
Ci > 2.6 being typical for ea rly-type galaxies and a < 2.6 being 
typical for late-type galaxies dStrateva et al. 2001 ). Our colour-cut 
has eliminated the disc galaxies with significant bulges, which is 
why we do not see a bias in the distribution of d across density 
groups. In addition, we still witness the steepening with density 
if we make a concentration index cut to select only galaxies with 
Ci < 2.6 (most galaxies in our sample have Ci < 2.6 to begin 
with). Since our sample primarily contains late disc type galaxies, 
we are not subject to bias from environmental effects related to the 
bulges of galaxies (there tend to be more bulges in disc galaxies 
in higher density environments). We would not detect in our study 
a trend in the TFR with environment which is really a tren d with 
morphological type, as in the results of Mas ters et alj d2006f) . 

• The environmental dependence of the TFR for late-type spi- 
rals is small (if any), meaning that the intrinsic properties of a 
galaxy play a principal role in driving the evolutio n of disc type 
galaxi e s. Our finding agrees with the earlier w ork by iMasters et alj 
d2006h . iMouhcine. Baldrv & Bamfor3 d2007l) found no environ- 
mental effect on the stellar mass-metallicity relation and conclude 
that the evolution of the galaxies is largely independent of their 
environments. They do find marginal increase in the chemical en- 
richment level at a fixed stellar mass in denser environments, so en- 
vironments can play a moderate role in galaxy evolution. An alter- 
ation in the TFR in denser environments could in theory be due to 
a temporary increase/accelerated star formation or to photometric 
or kinematic asymmetry due to dynamical interactions of galaxies. 
But we do not find strong evidence for this. 

• We have established a procedure for ascertaining if any trend 
in the SDSS-based TFRs with a variable of interest is due to physics 
or sample bias. In general, it is difficult to convert the observed 
SDSS-based TFRs as a function of dfwhm to TFRs as functions 
of Vcirc in a quick, simple manner because there is a large amount 
of scatter between wfwhm and v c i IC due to aperture biases, even 
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at higher redshifts where the apparent sizes of galaxies are smaller. 
Instead, it is most useful to work in magnitude versus Wpwhm space 
and convert an assumed TFR as a function of v c { rc into this regime 
with the aid of simulations of ti c i rc versus upwhm for the given size 
distribution of the galaxy sample. Looking at the SDSS line width- 
based TFR in multiple bands proves to be very useful because given 
two groups of galaxies which may have disparities in their v c i rc 
versus ufwhm relationship due to aperture biases, these disparities 
(e.g. steepening) should be reflected in the same way in the SDSS- 
based TFRs in all bands assuming the TFR is the same for both 
groups. 

• Using a steep er TFR for late type galaxies, like that of 
iMasters et all 12006). is also consistent with the data. The observed 
SDSS fFWHM-based slopes are still expected to steepen approxi- 
mately 7 per cent from the lowest to the highest density groups due 
to aperture bias effects. Two different sloped TFRs (as function of 
^circ), if their intercepts are carefully chosen, can predict the same 
slope and intercept for the observed SDSS VFWHM-based TFR. 
But they would also predict a statistically different distribution of 
vfwhm, w hich we do not obser ve ac ross density g r oups, so both 
the TFRs of IMasters et al.l J2006l) and lPizagno et alj J2007t) cannot 
be consistent at the same time across different density groups. The 
choice of a particular TFR does not affect our ability probe for dif- 
ferences in the inferred TFRs across density. 
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Figure Al. Coordinate system used for our model of a disc galaxy. The 
coordinates (r, <f>) describe a point in the plane of the disc and i is the ob- 
served inclination angle. The circle on the right hand side, perpendicular to 
the line-of-sight, is meant to indicate the 3" aperture. Coordinates (r', if)') 
are in the plane of the fibre, with the centre of the fibre (r' = 0) pointing to 
the centre of the galaxy (r = 0). 
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APPENDIX A: MODEL FOR SDSS-OBSERVED Ha 
PROFILE 

We describe here the simple model used to simulate synthetic popu- 
lations and understand the correlation between the SDSS-measured 
line widths and true disc velocities. In our model, we assume an 
infinitely thin disc with an inclination angle i to the observer's line- 
of-sight. The polar coordinates (r, ft) locate positions measured in 
the disc. A diagram of the coordinate system is presented in Fig- 
ure [Al] 

RCs of disc galaxies ca n often be well-modelled by an arctan 
function dPizagno et alj|2007h (although we consider the effect of 
rising and falling profiles as well in Appendix IB2t . The rotational 
velocity as a function from the centre of the disc is: 



2 ( r 

v(r) = — w c irc arctan — 
7r \rt 



(Al) 



where v c [ IC is the asymptotic circular velocity and r t is the turn- 
over radius (the distance where v(r) = u c i rc /2). Then, the ob- 
served velocity at (r, ft) is: 



VobsO", ft) = v(r) sin(i) cos(0) 



(A2) 



assuming we have subtracted off the systemic velocity of the 
galaxy. 

We assume a simple exponential surface brightness profile for 
the disc given by: 



I(r) = 7 exp 



(A3) 



where Iq is the central surface brightness and h is the exponential 
scale radius. The intensity observed from an infinitesimal element 
r dr dft is I(r)r dr d<f>. 

The spectrum observed from an infinitesimal element r dr d<f> 
will be centered at velocity w bs(r, (ft). We assume a Gaussian 
profile for the emission with velocity dispersion (standard devi- 
ation) a = 20 km s~ ( based on the typical value as obse rved 
for spiral type gala xies in lEpinat. Amram & Mar celin (20081) an d 
lEpinat et al.l ( l2010h . which may range between 10-30 km s _1 ). 
The dispersion a is the galaxy's velocity dispersion (for Ha) at 
an infinitesimal area in the disc, which we assume is constant as a 



function of r. The shape of the emission profile for element r dr < 
is then: 



P'(v) = I(r)r dr dcj> x N(v; v obB (r, (ft), a) 



(A4) 



where the notation Nix; fx, a) means a normal distribution as a 
function of x with mean fi and standard deviation a. 

To calculate the full Ha emission profile of the galaxy, the 
following integration can be performed, numerically, over the entire 
disc: 



Pfuii(v) = 



I(r)N(v;v obs (r,(p),a)r dr d(f>. (A5) 



To calculate the profile as observed by a 3" aperture, a similar 
integral can be evaluated numerically: 



Pfuii(v) = / / I{r)N(v;v oha (r,a>),a)F(r,ft)rdrd(f) 

(A6) 

where the added function F(r, (ft) describes the fraction of the 
emission from (r, (ft) that enters the 3" fibre. Establish new coor- 
dinates (r' t ft) to represent points in the plane P perpendicular to 
the line-of-sight with r — projecting onto r' — 0. Then the fi- 
bre can then be described as a circle C of radius 1.5" in plane P. 
Supposing that the fibre is centered on the galaxy's centre, r — 0, 
the fibre's centre is then at r' — 0. Without accounting for seeing, 
only points in the galaxy's disc that project to within circle C will 
contribute to the observed emission profile. A point (r, (ft) projects 
to plane P to have r' = proj(r, <ft) where 



proj(r, (ft) = r\J cos 2 {(ft) + sin 2 (</>) cos 2 (i). 

Thus, without accounting for seeing, the function F(r, q 

1 proj(r,0) < 1.5", 
proj(r,0)>1.5". 



Fxxjo seeing (^"7 ( 



(A7) 



(A8) 



To account for seeing, we assume a Gaussian convolution ker- 
nel withFWHM of 1.4". That is, the light that projects onto (r', ft') 
is redistributed as a 2D Gaussian. The fraction of this light that will 
enter the fibre is then the volume under the 2D Gaussian integrated 
over the surface with boundary C. This fraction only depends on 
r', the distance from the centre of the fibre, and not (ft . The integral 
is preformed numerically and the resulting fraction is shown in Fig- 
ure |A2] The fraction (blue line in Figure lA2l can be approximated 
as: 

f( r ') = - (Erf(1.5243 - 1.1157r') + Erf(1.5243 + 1.1157r')) 

(A9) 

where 



Erf(z) = -?=/ exp(-t 2 )dt. 

V7T Jo 



(A10) 



Hence, the function F(r, (ft), which takes into account seeing, 
is the following: 



F{r,ft) = /(proj(r, (ft)), 



(All) 
is given by Equa- 



where / is given by Equation |A9] and proj(r, 
tionfAT] 

Adding 1.4" seeing does not have a strong effect on the mea- 
sured FWHM of the emission profile. The difference between the 
FWHM with and without seeing is typically less than 2 — 5 km s" 1 , 
due to the fibre (3") being much larger than the seeing. This is use- 
ful to note because we have assumed a constant seeing while seeing 
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Figure A2. Figure showing the fraction of projected light from the galaxy 
to a distance r' from the fibre's centre that enters the 3" diameter fibre, 
without accounting for seeing (green) and with 1.4" FWHM seeing (blue). 

will vary in the SDSS but it does not have a big effect on the mea- 
sured line widths. 

An example Ha emission profile with and without the SDSS 
3" aperture bias is presented in Figure [T]for parameter values i = 
80°, h — 3", r t = 2 and v c i TC — 200 km s~ . The parameter 7o is 
just a normalization constant, which does not effect the calculation 
of the FWHM of the observed profile. 



comparable to varying the seeing (the change in the width is 
< 10 km s _1 ). And the variation in the observed width is virtually 
none for median-sized galaxies in our sample. This means that 
the SDSS fibre spectra are hardly sensitive to the shape of the RC 
beyond 0.8i? op t, where the profile may turn-over, remain flat, or 
rise. 

B3 Effect of Malmquist-type biases 

Malmquist-type biases may be present in any TFR sample. These 
biases may refer to the fact that dimmer galaxies at higher redshift 
could be missing from the sample due to the flux limit of the sam- 
ple, or flux errors scattering galaxies about the limit. These effects 
may we aken the observed sl ope of the TFR. The Tully-Fisher rela- 
tions of lPizagno et alj J2007I) suffers from Malmquist bias, leading 
them to underestimate the true slope of the TFR. Since we have not 
corrected for this here, our vfwhm TFR will also have a shallower 
slope, but, since our method is self consistent across the different 
environments, this weakening of the slope will not affect our con- 
clusions. 

This paper has been typeset from a TpX/ ETpX file prepared by the 
author. 



APPENDIX B: EFFECT OF SYSTEM ATICS ON 
SIMULATIONS 

Bl Effect of seeing on our results 

We have used a single seeing of 1.4" for all galaxies in our simu- 
lated of the observable fibre Ha line width. However, the spectro- 
scopic seeing for SDSS observations may range typically from 1" 
to 2". We run additional simulations of synthetic observations of 
25,000 galaxies with 1.1", 1.4", 1.7", 2.0" seeing. The relation 
between u c i rc and wfwhm also remains linear for galaxies with 
z > 0.045 assuming any of these seeing, and we find that in the 
worst-case, the slope of the relation between v c i IC and wfwhm is 
changed by 1.3 per cent from 1.1" to 2.0". 



B2 Effect of rotation curve shape on our results 

We have only considered the case where the RCs are asymptot- 
ically flat. However, disc galaxie s have been observed to have 
rising, flat, or falling RCs ( Catinella. Giovanelli & Havnesl 



rising , Hat, or falling Kl s (Latinella, Lriovanelli <x Haynes 
200rj- iPersic. Salucci & Stell Il996bl see also erratum 



Persic. Salucci & Stel 1996a). The RC profiles have been found to 



depend strongly on luminosity. Typically, the RCs may increase 
up to 0.8-Ropt (corresponding to 2.56ft for a pure disc profile) and 
then decline for the most luminous galaxies or continue increasing 
for the less luminous ones. Since we cannot reconstruct RCs from 
SDSS data, we consider the effects different RCs could have on 
the observed fibre widths in our sample in extreme cases and 
show that they are negligible. We consider a galaxy with apparent 
turn-over radius r t = 0.5" and h — 1.5", which is in the small 
end of our sample. We suppose that the RC has an arctan profile 
out to 0.8iiopt (2.56/i), followed by either a flat profile or an 
upward/downward profile with slope Av/Ar — 10 km s -1 ". 
But the variation in the observed velocity width is found to be 
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